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Optical excitations in a Bose gas are demonstrated to be very sensitive to many-body effects. At low 
temperature the momentum relaxation is provided by momentum exchange collisions, rather than 
by elastic collisions. A collective excitation mode forms, which in a Boltzmann gas is manifest in a 
collision shift and dramatic narrowing of spectral lines. In the BEC state, each spectral line splits 
into two components. The doubling of the optical excitations results from the physics analogous 
to that of the second sound. We present a theory of the line doubling, and calculate the oscillator 
strengths and linewidth. 



PACS: 

Recent experiments on the Bose-Einstein condensation 
in magnetically trapped gases @,|,Ul| stimulated stud- 
ies of many-body phenomena in such systems. Collec- 
tive excitations of several kinds Q have been investi- 
gated, including density waves [jl3| and the second sound 
0, excited by shaking the trapped BEC samples, or by 
an abrupt change of the interaction between particles 
jl4j]. All these responses are of hydrodynamical char- 
acter, with the excitation frequencies being not much 
greater than the trap frequency, and the wavelengths of 
order of the sample dimension. 

Optical excitations are also sensitive to the nature of 
quasiparticles and interactions between them. However, 
quasiparticles in a weakly non-ideal Bose gas, accord- 
ing to Bogoliubov theory, differ from free particles only 
at very low temperatures ||: ~ Tb = (AttTi 2 /m)aN <C 
3beC) where N is the gas density, and a is the s-wave 
scattering amplitude. In this letter we describe the 
many-body phenomena resulting from a change of the 
particle internal state due to optical transition. We shall 
see that these many-body effects can be significant even 
at temperatures as high as T ~ Tbecs such that the dif- 
ference of the quasiparticles and particles plays no role. 

The reason for the optical excitations being special is 
that the excited particle spends a long time in the state 
which is a coherent superposition of the ground state 
and the excited state. This time, which is of order of the 
inverse spectral line width, can be very long, so that the 
excited particle, while being in the superposition state, 
interacts with many other particles. The collisions of 
particles in different but non- orthogonal internal states 
have special properties, first discussed in the context of 
the problem of spin excitations in gases ||,|9|,|ll| (see also 
the work on collision frequency shifts in cesium fountains 

Let us consider a gas of identical atoms, all in the 
ground internal state |1). Suppose that one of the atoms 
is transferred to a superposition state |1'), having finite 



overlap with the ground state: |1') = a|2) + Gen- 
erally, there are three different kinds of collisions be- 
tween this atom and other atoms: (i) ordinary elastic 
collisions, occurring with the frequency r^ 1 = vtcl\ 2 N , 
where vt — (2T/m) 1 ^ 2 is thermal velocity, and ayi is 
the scattering length; (ii) forward collisions, in which the 
atoms interact without changing direction and velocity 
of motion; (iii) momentum exchange collisions, in which 
interacting atoms switch their trajectories and velocities. 
The forward and momentum exchange scattering rates 
are both given by T e ~* h = Airha^N/m, i.e. are temper- 
ature independent. Thus, at temperatures T < h 2 /ma 2 2 
the collisions (ii) and (iii) dominate. 

In this regime, named quantum gas by Bashkin B, col- 
lective effects appear due to collisions of type (iii). Note 
that forward collisions (ii) affect the non-interacting par- 
ticle picture only to the extent that they correct external 
potential seen by a particle by an amount proportional 
to the density of other particles. The role of the colli- 
sions (iii) is quite different. They provide mechanism of 
momentum transfer between different particles with the 
rate T o ~* h much higher than that of energy relaxation. 
Thus, at times < t c \ the particles can exchange momenta 
many times without changing the overall energy distri- 
bution. As a result, all excitations in the system with 
wavelength greater than urTexchj the path between mo- 
mentum exchange collisions, become collective modes. 

The collective phenomena in optical excitations arising 
due to momentum exchange collisions differ from those 
in the spin-wave dynamics [pp|JlO[|. The reason is that 
the interactions between atoms in different orbital states 
depend on the states, while for the atoms in different spin 
states all the interactions are the same. 
Problem: We consider optical excitations in a Bose gas of 
uniform density N . Each atom of the gas can be in the 
ground state or in the excited state (denoted by 1 and 
2, respectively). The interaction constants Xij, i,j = 
1,2 are related by the factor ATrh 2 /m with the s— wave 



1 



scattering lengths (an = a). 

The temperature range of interest, T < h 2 /ma 2 , in- 
cludes temperatures both above and below Tbec- For 
the collective effects to be pronounced, the excitation 
wavevector k should be < (uTT oxc h) _1 , which means that 
c\k\ <C wo, the excitation frequency. Thus, optimal situa- 
tion is that of a two-photon excitation with the photons' 
momenta equal and opposite. 

For example, this is realized in the recent study of the 
BEC in spin-polarized hydrogen by using the 15 — 2S 
Doppler free transition ]l5[ ]. In this experiment k ~ 
2ir/w, where w is the width of the laser beam used to 
excite the IS* — 2S transition. 

Below we study the excitation above and below Tbec- 
We calculate the frequency shift and the linewidth of 
the excitation. In the analysis we ignore spin degrees 
of freedom of the g well as the finite size of the gas 
sample confined in the trap. Thus, our results apply to 
the excitations probed near the central region of a shal- 
low trap, such that the density is nearly uniform. Also, 
our treatment is restricted to not very low temperatures, 
T > fi./T 0XC h ~ 2b, where one can ignore the difference 
between free particles and Bogoliubov quasiparticles. 
Summary of results: We show that the frequency of the 
excitation is shifted from that of a single atom by an 
amount different from the mean density shift. The reason 
for that is that the momentum exchange process couples 
different states degenerate in energy, and therefore, one 
has to apply degenerate perturbation theory to the many- 
body states to determine the correct frequency shift. For 
the Boltzmann gas, T > TbeCj by using the many-body 
version of the degenerate perturbation theory (the ran- 
dom phase approximation), we find the excitation fre- 
quency: 



wi^ 2 (fc) = w + 2(Ai2 - Xu)N + 
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3X 12 N 



(1) 



where the transferred momentum is assumed to be small 
in the sense discussed above: |fc| < X\ 2 N/vx- The res- 
onance frequency at k = varies linearly with density 
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. Such a behavior is observed in the experi- 
and the magnitude of the shift is found to be 
in agreement with (Q). 

The linewidth of the excitation (Q) due to momentum 
exchange collisions is found to be 



7exch ^ kv T e 



-(\ 12 N/kv T ) 2 



(2) 



(Elastic collisions contribute 7 c i — r el -C 7cxch-) The 
abrupt narrowing (B) of the spectral line at high densities 



N > kv T /X 
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(3) 



is similar to the motional (Dicke) narrowing effect (see 
Fig.|l|) This is because each momentum exchange colli- 
sion fully randomizes velocities, and so in the range (0) 



the collisions path, vx/Xi 2 N, becomes smaller than the 
wavelength 2i: jk. Another way to put it is to draw an 
analogy with Mossbauer effect, i.e., to say that in the 
regime (||) the recoil momentum is transferred collec- 
tively to all of the atoms interacting with the excited 
atom while it is in phase with the excitation field. 

Below Tbec, t ne main change in the spectrum is that 
the spectral line splits into two components, both shifted 
from the single particle frequency ujq. At k — the fre- 
quencies are given by 

wi.2 = w + (A12 — 2An)iV + Xi t 2 , 

(Xi, 2 - X n N c )(X lt2 - X 12 N T ) = X 2 12 N C N T , (4) 

where N c is the condensate density, and Nx = N — N c 
is the thermal density. The origin of the two frequencies 
(^) is that in the presence of the condensate the excited 
atom can be in two distinct states, "correlated" and "un- 
corrected" . In the correlated state the kinetic energy of 
excited atom is ~ T, and so it participates in momentum 
exchange collisions which shift w as discussed above. In 
the uncorrelated state, the velocity of the excited atom is 
zero, it does not take part in momentum exchange pro- 
cesses, and thus does not acquire the additional energy 
shift. Besides that, there are momentum exchange col- 
lisions which transfer the atom from the correlated to 
the uncorrelated state and backwards (see discussion be- 
low). The difference of the two frequencies (jj) gives the 
frequency of coherent oscillations between the two states. 

The behavior of the dispersion and the linewidth due 
to finite k as a function of temperature and density is 
qualitatively similar to that above Tbec (see Figs.|]||). 
Calculation The system is described by the Hamiltonian 
H = Hi + H 2 + Ti-ext- 
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Hi = E h a P a p + E 



An a pi a P2 a p3 a P4 



Pl+P3=P2 +P4 



n 2 = > \&- 



2m 



- w 



Pl+P3=P2+P4 



Ai2a+a P2 6+6 P4 



n 



bp +k a p + h.c. 



2m 

v 

E 

p 

p,k 

where a p and b p are canonical Bose operators for two 
internal states I and 2, w is the laser frequency, and A k 
are Fourier components of the laser field. We ignore the 
interaction A22 of atoms in the excited state, since we are 
interested in the absorption at small laser power, when 
the density of excited atoms is negligible with respect to 
total density N. 

As discussed above, the interesting regime is that of a 
(nearly) Doppler free two photon transition, k = ki + k 2 , 

w = wi + W2, where w ~ wo and |fc| ~ X\ 2 N/vt <C loq/c. 

The absorption spectrum is expressed through the two- 
particle correlation function by Kubo formula: 



I(w) = Im / e 



C(i),C+(0) 



dt 



(5) 
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(6) 



/j . k 



To evaluate the commutator in (g), we consider time 
evolution of the operators: idtrj Pt k = [%>,fc,7^i + H2]. In 
the situation of interest, when t c \ 3> T oxc h, the commu- 
tator can be evaluated by the RPA procedure |17[| . The 
resulting equations for rj p ,k(t) are linear: 



id t rj P ,k = ^ P ,kVp,k + Ai 2 w(p) ^2 rjp' 



(7) 



where 



(p + k) 



-^- + (Ai2-2A u )7V (8) 



To remind the reader how the RPA method works, 
consider the commutator 



V P ,k, ^2 a ti a P2 b P3 b : 

Pl+P3=P2+P4 

E 

Pl+P+k=p2+p. 



(9) 



a p a pi a P2 b pi 



E 



'p 3 °P4°p+k 



Pl+P3=P+Pi 



In the first term, according to the RPA procedure, we 
keep only two kinds of contributions: with pi = Pi, 
and with pi = p. Correspondingly, this gives Nrj Pi k + 
n (p)J2 P >Vp',k, which is the part of (0) proportional to 
A12. The second term in the commutator (|TJ) is of or- 
der of the inverse volume of the system, and so it can 
be neglected. (The first and second terms in (g) are of 
different magnitudes because there is only one particle 2 
in the systems, while particles 1 have finite density.) 

Solving the system ([?]), and plugging the result in (||), 
one gets the spectral power 



T(w) =ImJ2\A k \ 2 
fc 



i-Ai 2 n (fc,w) 



where 



n (A;,w) = 



n(p) 



d 3 p 



- a 



p ,k + i0 (27rft) 3 



(10) 



(11) 



The integrals in ([11]) over components of p normal to fc 
are evaluated explicitly, and the integral over longitudinal 
component of p is evaluated as a Cauchy integral, yielding 
an infinite (but quickly converging) sum over poles in the 
complex p plane p8| . The resulting spectrum is shown 
111 Fig.g for several increasing densities. 

To get the dispersion relation (|TJ) one expands IIo in 
small fc: kvx/XuN <C 1. The resulting dispersion (Q), 
which is much stronger than that for a free atom, is a 
collective mode effect ||,[l(| . 

To generalize the above calculation to the BEC situ- 
ation, we introduce ^b,fe — bi-ipo, where ^0 is the con- 
densate amplitude. In the RPA calculation the operator 



ffo,k is treated separately from other fj's. Now, the RPA 
dynamics is the following: 

id t r}p.k = Qp.kVp.k + Ai 2 n(p) rj . k + ^ Vp',k 

V V , 

id t rj ,k = (fi ,fc + AniV c )77o,fe + Ai 2 iV c ^ rjp'.k , (12) 

p' 

where N c = |"0o| 2 - Solving these equations yields the 
result of the form ([To|), 



(13) 



where 

U(k,Lu) = 



uj - fi , fc + (A12 - Xu)N c + iO 



(14) 



The spectrum at fc = consists of two sharp lines, with 
frequencies given by (0). 
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FIG. 1. Absorption spectrum QiaJ) , QXlJ) above BEC at fixed 
temperature and varying density: N = /-/Vbec, 0.5 < / < 1. 
The frequency shift Aui is defined relative to the that of a 
free atom at rest: Alu — uj — ujo — k 2 /2m. Spectral power X is 
normalized by particle density N. The excitation wavevector 
k is 0.5 in the units of Ai2.2Vbec /Vr- The interaction constant 
An is chosen to be 0. Note that the peak position follows the 
relation (Q). 

The physics of the doubling of the excitation below 
Tbec i s similar to that of the second sound. The excited 
particle can be in two different dynamical states, either 
having energy of order T and participating in momentum 
exchange processes, or being at rest. These two states 
are coupled by the transitions provided by momentum 
exchange processes (see (p"2|)), in which the excited par- 
ticle is transferred between the thermal and stationary 
states by exchanging with a condensate particle. The 
rate of such (coherent) transitions determines the fre- 
quency splitting in (H). 
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Relative strength of the spectral components given 
by ([l3|) depends on the condensate fraction N c /N. At 
N c -C N the peak due to the states with thermal ener- 
gies is much stronger than the other peak, which is due 
to atoms at rest. Because of that, we call the two peaks 
"normal" and "condensate", respectively. At iV c — * N 
the strengths of the two peaks become equal, and at 
T <C Ibec the peaks approach each other. This behavior 
is displayed in Figs.||,||. 




JU - Ato / (X N) 

FIG. 2. Absorption spectrum ©,(0) in the BEC regime, 
shown for temperature varying between and Tbec; density 
N fixed. The excitation wavevector k is 2/3 in the units of 
XhN/vt- Lines in the base plane indicate the peak positions 
(0) for k = 0. Note narrowing of the spectral line with de- 
creasing T, and strengthening of the condensate peak due to 
increasing condensate fraction. (The frequency shift Alj is 
defined in FigQ An = 0.) 




N /N, 



FIG. 3. Absorption spectrum in the BEC regime 

at fixed temperature and density varying from Nbec and up. 
Excitation wavevector k — 2.5 in the units of Ai2-/vbec/vt. 
Increasing condensate density leads to narrowing of the peaks 
and to strengthening of the condensate peak, as in Fig.jj]. ( Alj 
and T are defined in Fig.|l| An = 0.) 

Conclusion: We described phenomena in optical exci- 
tations resulting from cold collisions of the excited atom 
with other atoms of the gas. The main feature that makes 



the effect of collisions an interesting problem is that, 
during the optical transition, the excited atom spends a 
long time in the superposition of the ground and excited 
state. Because of many collisions that occur during this 
time, and because of their special coherent momentum- 
exchange character, the optical excitation represents a 
collective mode. The collective nature of the excitation 
is manifest in the frequency shift and in the linewidth 
narrowing. Below the BEC transition, due to the ex- 
change between the condensate and the normal compo- 
nent, the number of modes doubles, similar to the first 
and second acoustic modes. By using the RPA method, 
we calculate the temperature-dependent frequencies of 
the doublet and the oscillator strengths. 
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